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ABSTRACT 

The Phylogeny of Juniperus is presented based on nrDNA (ITS), petN-psbM, trnS-trnG, trnD-trnT, 
trnL-trnF sequencing (4411 bp) utilizing all currently recognized species. The major clades of the 
phylogenetic tree were congruent with the previously published phylogenetic tree of Mao et al. (2010) 
that had a subset of taxa of the current study. The lone species with serrate leaves in the eastern 
hemisphere, J. phoenicea, was found to be in a clade quite separated from the serrate junipers of North 
America. Juniperus phoenicea is referred to as 'pseudoserrate' to distinguish it from the semi-arid, serrate 
leaf junipers of the western hemisphere. Section Sabina is the most derived group and has radiated into 
niches in both the eastern and western hemispheres with approx. 60 species. 
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The genus Juniperus is comprised of approx. 75 species in 3 sections: 
sect. Caryocedrus , 1 species with large, blue, woody, 3-seeded cones, 
showing the fusing of 3 cone scales, with an Old World Mediterranean 
distribution (Adams, 2011, Adams and Schwarzbach, 2012a), 

sect. Juniperus, 14 species, 12 only in the eastern hemisphere, one endemic to 
North America and one species, J. communis , being circumboreal, seed cones 
blue or red, often with 3 seeds (Adams and Schwarzbach, 2012a) and 

sect. Sabina (approx. 60 species) with species about equally divided between 
the eastern and western hemispheres, seed cones with 1 to 13 seeds, blue, red- 
copper, rose, or brown (Adams, 2011). 

Section Sabina is divided into three major clades (Mao et al., 2010, 
Adams 2011): 

1. Serrate-leaf junipers of North America (21 species, Adams and 
Schwarzbach, 2011, 2013d), 

2. Turbinate-seed cones, single-seeded, entire-leaf junipers, eastern 
hemisphere (16 species, Adams and Schwarzbach, 2012b, 2013a, Zanoni 
and Adams, 1976, 1979) and 

3. Multi-seeded, entire-leaf junipers, both eastern and western hemispheres (23 
species, Adams and Schwarzbach, 2012c, 2013b). 
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The phylogenetic position of Juniperus (and Cupressaceae) in the plant kingdom (Fig. 1) shows 
Juniperus as a terminal clade and as one of the most advanced conifer genera (Rai et ah, 2008). Mao et 
al. (2010) demonstrated that the closest relatives of Juniperus are Cupressus (of eastern hemisphere) and 
the Hesperocyparis - Callitropsis - Xanthocyparis clade (Fig. 2). 


Although Mao et al. 

(2010) published a robust 
phylogeny of Juniperus, their 
principal purpose was to 
investigate the origins and 
evolutionary radiations of the 
major clades of Juniperus. As 
such, they utilized representative 
species from all clades, but not 
complete coverage of all known 
taxa. During the past few years, 
we have utilized DNA sequences 
to investigate the taxonomy of 
sections Caryocedrus and 
Juniperus (Adams and 
Schwarzbach, 2012a), section 
Sabina: serrate Juniperus of 
North America (Adams and 
Schwarzbach, 2006, 2012b, 

2013a, 2013c, 2013d; Adams and 
Nguyen, 2005); the turbinate 
seed cones, eastern hemisphere 
(Adams and Schwarzbach, 

2012b, 2013a) and the multi- 
seeded, entire-leaf junipers, both 
eastern and western hemispheres Figure 1. Maximum parsimony tree showing the position of 
(Adams and Schwarzbach, 2012c, Juniperus. Adapted from Rai et al. (2008). 

2013b). 

Bayesian tree (Mao et al. 2010) 
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Figure 2. Simplified Bayesian tree with genera collapsed showing Juniperus relationship to Cupressus and 
Hesperocyparis. Adapted from Mao et al., 2010. 


With the taxonomy and nomenclature having been addressed, the present report is to integrate these 
data into a robust phylogeny of Juniperus based on sequencing of nrDNA (ITS), petN-psbM, trnS-trnG, 
trnD-trnT and trnL-trnF including all known Juniperus species. 
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MATERIALS AND METHODS 

Specimens used in this study: see Adams and Schwarzbach (2012a, 2012b, 2012c, 2013a, 2013b, 
2013c). Voucher specimens are deposited at BAYLU herbarium Baylor University. 

One gram (fresh weight) of the foliage was placed in 20 g of activated silica gel and transported 
to the lab, thence stored at -20° C until the DNA was extracted. DNA was extracted from juniper leaves 
by use of a Qiagen mini-plant kit (Qiagen, Valencia, CA) as per manufacturer's instructions. 

Amplifications were performed in 30 pi reactions using 6 ng of genomic DNA, 1.5 units Epi¬ 
centre Fail-Safe Taq polymerase, 15 pi 2x buffer E (petN, trnD-T, trnL-F, trnS-G) or K (nrDNA) (final 
concentration: 50 mM KC1, 50 mM Tris-HCl (pH 8.3), 200 pM each dNTP, plus Epi-Centre proprietary 
enhancers with 1.5 - 3.5 mM MgCl 2 according to the buffer used) 1.8 pM each primer. See Adams, 
Bartel and Price (2009) for the ITS and petN-psbM primers utilized. The primers for trnD-trnT, trnL-trnF 
and trnS-trnG regions have been previously reported (Adams and Kauffmann, 2010). 

The PCR reaction was subjected to purification by agarose gel electrophoresis. In each case, the 
band was excised and purified using a Qiagen QIAquick gel extraction kit (Qiagen, Valencia, CA). The 
gel purified DNA band with the appropriate sequencing primer was sent to McLab Inc. (San Francisco) 
for sequencing. Sequences for both strands were edited and a consensus sequence was produced using 
Chromas, version 2.31 (Technelysium Pty Ltd.) or Sequencher v. 5 (genecodes.com). Sequence datasets 
were analyzed using Geneious v. R6-1 (Biomatters. Available from http://www.geneious.com/ ), the 
MAFFT alignment program. Further analyses utilized the Bayesian analysis software Mr. Bayes v.3.1 
(Ronquist and Huelsenbeck 2003). For phylogenetic analyses, appropriate nucleotide substitution models 
were selected using Modeltest v3.7 (Posada and Crandall 1998) and Akaike's information criterion. 

RESULTS AND DISCUSSION 

The phylogeny of Juniperus utilized the most informative gene regions (nrDNA, petN-psbM, 
trnSG, trnDT and tmLF). The tree is similar (Fig. 3) to Mao et al. (2010), except the positions of J. 
californica , J. durangensis, and for their clade IV which is now resolved into excelsa and chinensis 
groups (Fig. 3). The use of duplicate samples for most taxa appears to have stabilized the Bayesian tree 
in many places increasing branch support. In addition, all known taxa are included (approx. 100), 
compared with 51 taxa by Mao et al. (2010) and this has strengthened the tree. 

Several 'problem' taxa present difficulties: J. phoenicea and J. p. var. turbinata stand loosely 
affiliated with sect. Sabina (Fig. 3). These taxa have small serrations on the leaf margins, but are denoted 
as 'pseudoserrate' (Fig. 3). It seems unlikely that serrate leaf margins in the eastern and western 
hemispheres is a homologous character, but has arisen independently as J. phoenicea is not in the clade 
with the serrate, semi-arid junipers of the western hemisphere. Also, J. erectopatens and J. microsperma 
form an unusual clade that does not nest into the J. chinensis clade (Fig. 3). Juniperus ashei and J. a. var. 
ovata, now J. ovata (R. P. Adams) R. P. Adams, are in separate clades (see Adams, 2008; Adams and 
Schwarzbach, 2013d for discussion). 

Mao et al. (2010) used three Juniperus fossil dates: J. pauli (ca. > 33.0 mya, cf. extant J. sabina 
and allies), J. creedensis (ca. > 23.0 mya, cf. J. californica / J. osteosperma ), and J. desatoyana (ca. > 
16.0 mya, cf. J. occidentals / J. osteosperma ). They postulated the serrate, semi-arid junipers migrated 
from the eastern to the western hemisphere via the North American Land Bridge (NALB) ca. 47 - 30.3 
mya (Fig. 4). 
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Figure 3. Bayesian tree of the Phylogeny of Juniperus. Numbers on lines are posterior probabilities. 




















































































Phytologia (May 2013) 95(2) 


183 



Figure 4. Possible migration of ancestor of the serrate, semi-arid junipers from the Mediterranean to the 
western hemisphere via the NALB. 


The fossil J. creedensis of the Creede geoflora 
(ca. > 23.0 my a) bears a striking resemblance to 
present-day J. californica (Fig. 5). Because the 
present-day J. californica appears little changed from 
the fossil, J. creedensis , it may be that the serrate 
junipers in North America are much older than 
thought. It might be noted that Axelrod (1987) 
described a second juniper from the Creede geoflora 
as J. gracillensis that he thought was similar to extant 
J. flaccida, but Wolfe and Schorn (1990) have 
identified the specimen as Eleopoldia lipmanii 
(Rosaceae). 



J. creedensis 


J. californica 


Fig. 5. Juniperus creedensis Axelrod paratype 
and present-day J. californica. 


The Madrean-Tethyan vegetation belts in Eurasia and North America may have been continuous 
during the Eocene and Oligocene (Axelrod, 1975, Wen and Ickert-Bond, 2009), such that Juniperus 
section Sabina might have had a wider distribution (Fig. 6). So it is possible that the serrate-leaf junipers 
may have existed in the Madrean-Tethyan vegetation belts in both Eurasia and North America during the 
same period (Fig. 6), and there may have been exchanges via the North Atlantic Land Bridge (NALB). 
However, one should note the Madrean-Tethyan vegetation, depicted in Figure 6, predates the ages of any 
known juniper fossils. 



Figure 6. Madrean - Tethyan vegetation zones of Axelrod (adapted from Axelrod (1975) and Wen and 
Ickert-Bond (2009). 
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It is unfortunate that the only serrate-leaf species (. J . phoenicea) that is extant in the eastern 
hemisphere, has DNA so different that it is poorly grouped with any clade (Fig. 3). It is quite removed 
from the serrate junipers clade (Fig. 3). At present, it seems appropriate to consider J. phoenicea as 
'pseudoserrate' and of a different lineage than the serrate junipers of North America. If the serrate leaves 
of J. phoenicea are not homologous to the serrate leaves of junipers in North America, then we are left 
with no extant (or known fossils) of truly serrate junipers in the eastern hemisphere. 

The migration of the smooth leaf members of sect. Sabina to the western hemisphere is thought to 
be more recent (17.6-5.5 mya, Mao et al., 2010) and those dates are younger than the fossil J. creedensis 
of the Creede geoflora (Axelrod, 1987). Since J. phoenicea does not appear to be a true member of the 
serrate junipers and no serrate juniper fossils have been found in the eastern hemisphere, the serrate 
junipers may be endemic to the western hemisphere. Undoubtedly, additional fossils will be found some 
day to help resolve the question. 

Mao et al. (2009) argues that the movement of sect. Sabina to the western (17.6-5.5 mya) is too 
young for migration across the North Atlantic Land Bridge (NALB), but possible via the Bering Land 
Bridge (BLB). Because sect. Sabina species such as J. sabina and J. davurica are quite cold adapted, 
they could have migrated to produce the ancestors that gave rise to the current, cold climate, western 
hemisphere species such as J. horizontalis and J. scopulorum. Juniperus davurica is the northeastern- 
most species in northeast Asia (in sect. Sabina ) and could have provided ancestral stock to migrate across 
the Bering Land Bridge (Fig. 7). Notice that J. davurica - J. sabina are in a sister clade to the smooth- 
leaf juniper of North America (Fig. 3), supporting the concept of migration from northeastern Asia via the 
BLB (Fig. 7.) 



Figure 7. Possible migration pathway of the smooth leaf junipers of North America. 


Juniperus communis is an interesting taxon in that it is the most weedy (invasive) species in sect. 
Juniperus. Its seed cones are especially juicy and attractive to birds. It is found in disturbed habitats as 
an invasive weed in Hungary and central Europe as well as in North America (as J. c. var. depress a). 
Juniperus communis and J. c. var. depressa form a boreal distribution in the higher latitudes around the 
northern hemispheres. Adams and Schwarzbach (2012a) examined the taxonomy of J. communis and 
found it to be very complex. The species is comprised of several morphological varieties that are closely 
linked by only a few mutations (Fig. 8). Notice the Kamchatka group (Fig. 8) is closely linked (6 
mutational events, MEs) to J. communis var. nipponica, Japan, thence to J. c. var. megistocarpa (NA, 5 
MEs). 


The North America communis group is equally linked between the Japan and Europe-Central Asia 
groups. Thus, the linkage map gives equal support to the Bering Land Bridge and North Atlantic island 
hopping model for the origin of J. communis in North America. The situation was previously more 











Phytologia (May 2013) 95(2) 


185 


unclear when J. jackii was included in J. communis ( J . c. var. jackii). However, J. jackii is clearly quite 
differentiated (20 MEs from J. mairei, Gansu, China; 21 MEs from J. c. var. megistocarpa, NA and Fig. 
3), but the data is equivocal as to whether its origin is from the BLB or North Atlantic island hopping 
model. It grows on serpentine and volcanic basalt of quite recent origin in the Cascade Range of western 
Oregon/ northern California. 


Minimum Spanning Network 
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Blue Seed Cone Juniperus, 
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Figure 8. Minimum spanning network of blue seed cone junipers. Numbers on the links are MEs 
(mutational events). Adapted from Adams and Schwarzbach (2012a). 


A diagrammatic representation of the possible migrations of J. communis (and J. jackii) is shown 
in Figure 9. The migration dates proposed by Mao et al. (2010) seem consistent with the recent habitat 
availability for J. jackii and support the observed lack of differentiation among morphological varieties of 
J. communis (Adams and Schwarzbach, 2012a). 



Figure 9. Possible migrations of J. communis to North America. Based on data in Mao et al. (2010). 
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Wen and Ickert-Bond (2009) summarized data from 17 studies concerning Madrean-Tethyan 
disjunctions. Their summaries are useful in the present discussion. They concluded (Fig. 10) that: 53% 
of the inter-continental migrations was by the North Atlantic Land Bridge; 40% was by long distance 
dispersal and 7% by the Bering Land Bridge (BLB). Their summary of the directional data indicated the 
origins as: 86% from eastern to western hemisphere; 7% from western to eastern hemisphere and for 7% 
the direction was uncertain (Fig. 10). This trend broadly supports the conclusions of Mao et al. (2010) 
and the present study. 



Summary of migration patterns 
(Wen & Ickert-Bond, 2009) 



direction 
uncertain, 7% 

w. hemisphere 
to e. hemisphere 
7% 


e. hemisphere to 
w. hemisphere 
86 % 


Summary of direction of origin 
(Wen & Ickert-Bond, 2009) 


Figure 10. Summary of trends based on 17 studies of Madrean-Tethyan disjunctions. Adapted from Wen 
and Ickert-Bond (2009). 


CONCLUSIONS 

The present phylogenetic results are compatible with the results of Mao et al. (2010). In the 
present report, the clades are better defined and clade IV (Mao et al. 2010) is now resolved into three 
clades of the J. excelsa group and the J. chinensis group, with J. erectopatens and J. microsperma 
forming a somewhat intermediate clade. The lone species with serrate leaves in the eastern hemisphere, 
J. phoenicea, was found to be in a clade quite separated from the serrate junipers of North America. It 
appears that the evolution of serrate leaves occurred independently in the eastern hemisphere. Juniperus 
phoenicea is referred to as 'pseudoserrate' to distinguish it from the semi-arid, serrate leaf junipers of the 
western hemisphere. Section Sabina is the most advanced group and has radiated into niches in both the 
eastern and western hemispheres with approx. 60 species. Additional fossils are needed from older 
formations to clarify the evolution of the genus. 
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